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As a sister work of Ref . [3| , we incorporate the color-screening effect due to light quark pair creation 
into the heavy quark-antiquark potential, and investigate the effects of screened potential on the 
spectrum of bottomonium. We calculate the masses, electromagnetic decays, and El transitions 
of bottomonium states. We find that the fine splittings between Xb-J (J=0,l,2) states are in good 
agreement with experimental data, and the El transition rates of T(25') ^ IXbji^P) and T(35') 
7X6j(2P) (J=0,l,2) all agree with data within experimental errors. In particular, the mass of T(6S') 
is lowered down to match that of the T(11020), which is smaller than the predictions of the linear 
potential models by more than 100 MeV. Comparison between charmonium and bottomonium in 
some related problems is also discussed. 
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I. INTRODUCTION 

Potential models have been successful in describing the 
spectra below the open-flavor thresholds for both char- 
monia and bottomonia. However, it is well-known that 
these quenched potential models, which incorporates a 
Coulomb term at short distances and a linear confining 
potential at large distances 0, [Ij 0] i '^i^l overestimate the 
masses of heavy quarkonia above the open-flavor thresh- 
olds. Some distinct examples are the Xc2'/^(3930) in 
the charmonium system and the T(6S') in the bottomo- 
nium system, of which the observed masses are about 50 
and 90 MeV respectively lower than that predicted by 
the typical relativized potential model of Godfrey and 
Isgur Q, and even 120 MeV lower for T{6S) than the 
prediction of the Cornell model 0, Q . 

This is probably because, the linear potential, which is 
expected to be dominant at large distances, is screened or 
softened by the vacuum polarization effect of the dynami- 
cal light quark pairs Q . This screening or string breaking 
effect has been demonstrated indirectly by the investiga- 
tion of the mixing of static heavy quark-antiquark (QQ) 
string with a static heavy-light meson-antimeson {Qq- 
Qq) system in the nj = 2 lattice QCD calculations p, 
and has also been implied recently by the calculations 
within some holographic QCD models 

However, since the simulations of lattice QCD still have 
large uncertainties and difficulties in handling higher ex- 
cited states, it should be useful to improve the potential 
model itself to incorporate the screening effect and com- 
pare the model predictions with experimental data, as a 
phenomenological way to investigate the screening effects 
on heavy quarkonium spectrum. 

Such screened potential models [1, d, [13, [ll| were pro- 
posed many years ago in the study of heavy quarkonium 
and heavy flavor mesons, as well as light hadronsfl^. 
The main feature of these screened potential models in 



the spectrum is that the masses of the higher excited 
states are lowered. 

In recent years the screened potential models have 
again been used to investigate the heavy quarkonium 
spectrum and leptonic decay widths [l^. In Ref.|l| we 
have reinvestigated the charmonium spectrum within the 
screened potential model suggested by Chao and Liu [1| 
and assigned some newly discovered charmonium-like res- 
onances as conventional higher charmonium states. 

On the experimental side, aside from abundant res- 
onances discovered recently in the charmonium region, 
progress in the bottomium region has also been made. 
The T(1Z?) was observed by CLEO collaboration \\M in 
2004 and the rji, was observed by BaBar collaboration [1^ 
in 2008. One may expect more bottomium states will be 
observed in the future by BaBar, Belle and CLEO. So 
it is important to reinvestigate the bottomonium system 
within the screened potential model. 

In this paper, as a sister work of we calculate the 
mass spectrum and electromagnetic decay and transition 
rates of bottomonium especially the higher bottomonium 
using a non-relativistic Schrodinger equation with the 
Coulomb potential plus a screened linear potential, which 
is nearly the same as that in The model predictions 
are similar to that of [l| for charmonium. The mass of 
T{QS) is lowered to be consistent with its experimental 
value. 

In the following, we first introduce the screened po- 
tential model in Sec. II, and then study the mass spec- 
trum and decay and transition processes of bottomonia 
in Sec. III. In Sec. IV we will discuss some features of our 
result for the bottomonium states. A summary will be 
given in Sec. V. 
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II. THE SCREENED POTENTIAL MODEL 



For the model parameters, we take 



As a minimal model describing the bottomonium spec- 
trum we use a non-relativistic potential model with the 
screening effect being considered as in We use a po- 
tential as 



Vscrir) ^Vvir) + Vsir), 



Vv{r) 



Vs{r) = A( 



4 ac 

1 - e-'"' , 



C. 



(1) 



(2) 



(3) 



Here /x is the screening factor which makes the long-range 
scalar potential of Vscrir) become flat when r ^ and 
still linearly rising when t ^ -j^, and A is the linear po- 
tential slope (the string tension), which is taken to be the 
same as for charmonium^. VV(r) represents the vector- 
like one-gluon exchange potential, ac is the coefficient 
of the Coulomb potential. C is a constant related to the 
normalization of energy levels of the QQ system. 

The spin-dependent interactions include three parts as 
follows. The spin-spin contact hyperfine interaction is 

Hss^^^l{r)S,-Sj,, (4) 



ac = 0.37, as = 0.18, 
fi = 0.056 GeV, C = 0.677 GeV, 
mb ^ 4.4 GeV, a ^ 3.3 GeV, X = 0.2106^^ (10) 

where ac ~ as{mbVb) and as ~ as{2'mb) are essentially 
the strong coupling constants but at different scales. The 
former is for large distances and used to determine the 
spectrum while the latter is for short-distances and used 
for QCD radiative corrections in bottomonium decays 
(see below in next section). 

Here ^ is the characteristic scale for color screening, 
and l//i is about 3.5 fm, implying that at distances 
larger than l//x the static color source in the bb sys- 
tem gradually becomes neutralized by the produced light 
quark pair, and string breaking emerges. Note that 
here fi is smaller than that of charmonium in [l|, where 
H = 0.0979 GeV^ corresponding to 2 fm. In Sec. IV we 
will discuss the reason for the difference in /i between bb 
and cc systems. 

With these values of the parameters for the potential, 
we can calculate the spectrum of bottomonium. The re- 
sults are shown in Table U For comparison, we also list 
the experimental values and those predicted by the 
linear potential model [3| in Table HI 



where Sa{r) is usually taken to be S{f) in nonrela- 
tivistic potential models, but here we take Sa-(r) = 
{a I ypIxY e^° as in Ref. [lB| since it is an artifact of 
an 0(w^/c^) expansion of the T-matrix [17] in a range 
comparable to 1/mf,. 

The spin-orbit term and the tensor term take the com- 
mon forms 

Hi^s = ^{Wy{r)-V's{T))L-S, (5) 



2mfr 



and 



l_(il/(r)-VV(r))r. 



12m^ r 



(6) 



These spin-dependent interactions are dealt with pertur- 
batively. They are diagonal in a \J,L,S > basis with the 
matrix elements 



< Sb- Sb >- -S'^ - -, 



(7) 



<L-S>^ [J{J + 1) - L{L S{S + l)]/2 (8) 

and the tensor operator T has nonvanishing diagonal ma- 
trix elements only between L > spin-triplet states, 
which are 



L 

'6(2L+3) 



1 



<^Lj\T\^Lj >= < 



1 

6 ' 



(L + 1) 
6(2L-1) ' 



J = L 
J = L 
J = L- 1 



(9) 



III. SOME DECAY PROCESSES 

A. Leptonic decays 

The electronic decay width of the vector meson is given 
by the Van Royen-Weisskopf formula Il9l with QCD ra- 
diative corrections taken into account [20|. 



TeeinS) 



n.(0)P(l~^^), (11) 



TeeinD) 



25a'el „ 



16 as 



), (12) 



where MnsiMno) is the mass for nS{nD), e;, = i is the 
6 quark charge, a is the fine structure constant, i?„s(0) 
is the radial S wave- function at the origin, and i?„£)(0) 
is the second derivative of the radial D wave-function at 
the origin. 

With the chosen parameters (fTU]) . we get the results 
that are listed in table |TT1 We also list other two groups' 
results [2l|, for comparison. We can see that our re- 
sults are consistent with the experimental data. 



B. Two-photon decays 

In the nonrelativistic limit, the two-photon decay 
widths of the ^5*0, ^Pq, and states can be written 
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as [2E 



r^Y5o-77) 
r^«(^P2->77) 



3a\t\Rnsi 



36a\t\Kpm' 
5m^ 



(13) 
(14) 
(15) 



The first-order QCD radiative corrections to the two- 
photon decay rates can be accounted for as [2^ 



TCSo 


77) 




77) [1 - 


as , 

TT ^ 3 






77) 


= rNR (3po - 


-4 77) [1 - 


TT 




r(3p2 


77) 


= rN« (3p2 - 


-4 77) [1 - 


16 as, 

' 3 TT ^' 


(18) 


We 


can see 


that T(}So - 


-> 77) oc 


|i?nS(0)P, 


which arc 



sensitive to the details of the potential near the origin. 
So we take 



T{% -> 77) ^^'^"i^JJ^ ririnS) 

1 ee(nS) 



(19) 



to eliminate this uncertainty. 

In the nonrelativistic limit, we can also replace mt by 
M/2, where M is the mass of the corresponding bottomo- 
nium state. The results are listed in Table Till Predic- 
tions of some other models (see Refs. d, |24 l25l. [26l , [27l . 
[2^, [HI) are listed for comparison. We can see our results 
are a bit larger than most models but are consistent with 

Refs. M,m- 



C. El transitions 

For the El transitions within the bottomonium system, 
we use the formula of Ref . [301 : 



,2S + lLj 



n 



/ 2S' + 1t / 



Lj'+7) 



rEi(n^ 

-^Cf^Sss' ela\{f\r\^) p 

where E^ is the emitted photon energy. 
The spatial matrix element 



< /|r|i >= / Rf{r)R,{r)r^dr , 
Jo 



(20) 



(21) 



involves the initial and final state radial wave functions, 
and the angular matrix element C/^ is 



Cf, = max(L, L') (2J' + 1) | j ^ 1 



V J'S 



(22) 



Our results are listed in Table IVIl The widths calcu- 
lated by the zeroth-order wave functions are marked by 



SNRq and those by the first-order relativistically cor- 
rected wave functions are marked by SNRi. 

For the first-order relativistic corrections to the 
wave functions, we include the spin-dependent part of 
(HI),®,® and the spin-independent part as [3l| 



1 



SI 



pi 



2m 



o 2 

2mg 



(23) 



where Pi and P2 are the momenta of b and b quarks 
in the rest frame of bottomonium, respectively, which 
satisfy Pi = — P2 = P, 9 is the unit second-order tensor, 
and {{ }} is the Gromes's notation 

{{A-n-B}} = ^{AB : n+A-m+B-^A+n : AB), (24) 

where 5R is any second-order tensor. 

Note that we do not include the contributions from 
the scalar potential in Hgj since it is still unclear how to 
deal with the spin-independent corrections arising from 
the scalar potential theoretically. 

We also fist in Table [Vlj the resuhs of Ref. 30] which 
uses a potential obtained from the inverse-scattering 
method for comparison. 

Both the SNRq and Rcf. results of El transitions 
are larger than most of the experimental values, but we 
see that in SNRi the predicted widths get decreased and 
fit the experimental values quite well as long as the first- 
order relativistic corrections to the wave functions are 
taken into account. 



IV. DISCUSSIONS 

A. T(11020) 

T (11020), the candidate of T(6S), was observed in 
e+e^ annihilation in 1985 [32,[33|- Its PDG mass and 
full width ;i8| are 



M = 11019±8AfeV, 
r = 79 ± 16 MeV. 



(25) 



Recently, the BaBar collaboration [3J] has remeasured 
the e+e^ — s- bb cross section by an energy scan in the 
range of 10.54 GeT^ to 11.20 Gel^ and get the mass and 
full width as 



M = 10996 ± 2 MeV, 
r = 37 ± 3 MeV. 



(26) 



Despite of the discrepancy in the mass and full width 
jiven by BaBar and PDG, the observed mass is much 
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smaller than that predicted by the linearpotential mod- 
els. For example, the Cornell model [3, Q predicted 
11.14 Gey, which is 121(144) MeF larger than the ex- 
perimental value of PDG(BaBar), and the modified Cor- 
nell model [33] gives 11.113 GeF, which is 104(127) MeV 
larger than the experimental value of PDG(BaBar). The 
relativized potential model of Godfrey and Isgur 4] gives 
11.10 Gey, which is 91(114) MeV larger than the exper- 
imental value of PDG(BaBar). 

Evidently, the mass of T(6S') is overestimated by the 
quenched potential models by more than 100 MeV. If we 
take the screening effect into account, we find, in our 
model, the mass of T(6S') to be 11.023 GeT^, which is 
very close to the observed value of PDG(BaBar). The 
consistence of our predicted mass with the experimental 
value of T(11020) indicates the significance of the screen- 
ing effect on higher excited bottomonium states. 



B. Hyperfine and Fine Splittings 

We use Eq. ([4]) to calculate the hyperfine splittings be- 
tween T(nS) and r]ii{nS), where S^ir) is taken to be 
Sa{r) = {a / y/ir)^ e^"' ^ as in Ref.[l6]. a has a mag- 
nitude of order mq and it represents some relativistic 
smearing effects [l6j. 

We choose the observed splitting between J/ip and rjc 
as input to determine a for charmonium and have ob- 
tained a good fit for the i/'(3686)-77c(3637) splitting [l[). 
Here we use the observed rjb P^] and T(IS') masses to de- 
termine a for bottomonium, and find a to be 3.3 GeV in 
our model [see pO))]. The hyperfine splittings for charmo- 
nium and bottomonium systems are listed in Tab|Vl In 
comparison, we also list the results of RefJSGj] , which uses 
the Buchmiiller-Tye potential and Ref.jJ], which uses a 
relativized funnel potential. 

We also list the results of the fine splittings, which are 
calculated by using ([5]) and (O, between P-wave multi- 
plets for both charmonium and bottomonium in the same 
table. We can see that our results fit the experimental 
values quite well and are compatible with Ref.[36j and 
Ref.[3. 



C. El transitions 

We have calculated the El transition widths for bot- 
tomonium using the zeroth-order wave functions, which 
are marked by SNRq, and the first-order relativistically 
corrected wave functions, which are marked by SNRi. 
The results are listed in Table IVII along with the results 
from the potential model in Rcf. 30] , in which the poten- 
tial is determined by the inverse-scattering method, for 
comparison. 

We find our results are compatible with experimen- 
tal values and Ref. [lO] for most channels. Relativistic 
corrections to the wave functions tend to reduce the 



El transition widths for most channels and give bet- 
ter fit with experimental values. Note that for the 
T(3S') -> 7Xho transition our result, 0.07(0.05) KeV with 
the zcroth(first)-order wave functions, is in agreement 
with the experimental value 0.061 ± 0.023 KeV, while 
that of Ref.[33| (0.007 XeF) is too smah. 

But our calculated transition widths for T(3S') 
lXbj{J — lj2) are too large as compared with exper- 
imental data. These may indicate that for the radi- 
ally suppressed El transition widths (e.g. T(3S') 
1Xbj{J = 0,1,2)) the theoretical values are very sen- 
sitive to model details, and further improvement for the 
model and the calculation is needed. 



D. Screening parameter fj, 

We find the screening parameter /z, which represents 
the energy scale related to the creation of the Qq and 
Qq pair or the distance when that beyond r ~ 1/^ 
the screening effect becomes important, is smaller for 
bottomonium (/i = 0.056 GeV) than that for charmo- 
nimn (p = 0.0979 Gey)[l], if we try to fit the bottomo- 
nium spectrum. We need to understand this difference 
of fi between the bb and cc systems. 

It is known that the string breaking is due to the cre- 
ation of light quark pairs, i.e., the formation of Qq-Qq 
mesons. Note that for the bb system the energy differ- 
ence between the bb bound state T(IS') and the open 
bottom threshold of BB meson pair is 1.1 GeV, whereas 
for the cc system the energy difference between the cc 
bound state J/ip and the open charm threshold of DD 
meson pair is only 0.63 GeV. This implies that for the bb 
system more energy needs to be stored (or equivalently 
a longer flux tube is neded) before the string breaking 
occurs than that for the cc system. 



V. SUMMARY AND CONCLUSIONS 

In this paper, as a sister work of [l|, we incorporate 
the color-screening (string breaking) effect due to light 
quark pair creation into the heavy quark-antiquark long- 
range confinement potential, and investigate the effects 
of screened potential on the spectrum of bottomonium. 
We calculate the masses, electromagnetic decays, and El 
transitions of bottomonium states in the nonrelativistic 
screened potential model. 

We find that the screening parameter /i is smaller for 
bottomonium than that for charmonium if we try to fit 
the bottomonium spectrum, and this may be understood 
as due to the difference between bb and cc systems in the 
energy to be stored before the string breaking occurs. 
The masses predicted in the screened potential model 
are considerably lower for higher bottomonium states. 
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compared with the unscreened potential model. Espe- 
cially, the mass of T(6S') is lowered down to match that of 
T(11020), whereas the hnear potential model predictions 
are more than 100 MeV higher than the experimental 
value. The fine splittings of P-wavc bottomonium states, 
and El transition rates and leptonic decay widths are 
found to be compatible with experimental data within 
errors. 

We hope our investigation for the bottomonium system 
with screened potential model will be useful in the future 
study of bottomonium physics. 
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TABLE L Experimental and theoretical mass spectrum of bottomonium states. The experimental masses are PDG [3 averages. 
The masses are in units of MeV except for Ref. 4] which is in GeV. The averaged radiuses are in units of fm. The results of 
our screened potential model are shown in comparison with that of Ref. 3. 



State 


Expt. 


Theor. 


of ours 


Theor. of Ref [4] 






Mass 




Mass 


IS T(l^Si) 


9460.30 ± 0.26 


9460 


0.23 


9.46 




9388.9l2:3 ± 2.7 


9389 




9.40 


2S T(2^Si) 


10023.26 ± 0.31 


10016 


0.52 


10.00 






9987 




9.98 


3S T(3^Si) 


10355.2 ±0.5 


10351 


0.78 


10.35 






10330 




10.34 


4S T(4^Si) 


10579.4 ± 1.2 


10611 


1.02 


10.63 


r?5(4^So) 




10595 






5S T(5^Si) 


10865 ± 8 


10831 


1.24 


10.88 


?7b(5^So) 




10817 






6S T(6^Si) 


11019 ±8 


11023 


1.45 


11.10 


??b(6'So) 




11011 






7S T(7^Si) 




11193 


1.66 




J76(7^So) 




11183 






IP X62(1"P2) 


9912.21 ±0.26 ±0.31 


9918 


0.42 


9.90 


Xf>i(l'Pi) 


9892.78 ±0.26 ±0.31 


9897 




9.88 


X6o(l'Po) 
/li,(l'Pl) 


9859.44 ± 0.42 ± 0.31 


9865 




9.85 




9903 




9.88 


2P Xf.2(2^P2) 


10268.65 ± 0.22 ± 0.50 


10269 


0.69 


10.26 


Xw(2^Pi) 


10255.46 ± 0.22 ± 0.50 


10251 




10.25 


Xw(2^Po) 


10232.5 ±0.4 ±0.5 


10226 




10.23 


/ic(2^Pi) 




10256 




10.25 


3P X62(3^P2) 




10540 


0.93 




Xw(3'Pi) 




10524 






Xw(33Po) 




10502 










10529 






4P Xb2(4^P2) 




10767 


1.15 




Xw(4'Pi) 




10753 






Xw(43Po) 




10732 






/ii,(4^Pi) 




10757 






5P Xf>2(5"P2) 




10965 


1.37 




Xi>i(5^Pi) 




10951 






Xf.o(53Po) 




10933 








10955 






ID V3(1^D3) 




10156 


0.57 


10.16 


V'2(1^D2) 


10161 ±0.6 ± 1.6 


10151 




10.15 






10145 




10.14 


J7c2(l'D2) 




10152 




10.15 


2D V3(2^D3) 




10442 


0.82 


10.45 


^'2(2^02) 




10438 




10.45 


t/'(23Di) 




10432 




10.44 


r;c2(2^D2) 




10439 




10.45 


3D ^3(3^03) 




10680 


1.05 




V'2(3^D2) 




10676 






V'(3^Di) 




10670 






»7c2(3^D2) 




10677 






4D t/>3(4^D3) 




10886 


1.27 




^'2(4^02) 
t/'(3''Di) 




10882 








10877 






r;c2(4^D2) 




10883 






5D ^3(5^03) 




11069 


1.49 




V'2(5^D2) 




11065 






V'(5'Di) 




11060 






»7c2(5^D2) 




11066 
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TABLE II: Leptonic decay widths (in units of KeV) for bottomonium states in our screened potential model. The widths 
calculated with and without QCD corrections are marked by Fee and F^e respectively. The experimental values are taken from 
PDG [18jj. Predictions of two other models f2ll, [25| are listed for comparison. 



state 


pU 

^ ee 


Tee 


Ref.f21] 


Ref.f22] 


ExpflS] 


1'' Si (9460) 


2.31 


1.60 


1.314 


1.320 


1.340 ± 0.018 


2^ Si (10023) 


0.92 


0.64 


0.576 


0.628 


0.612 ± 0.011 


3^ Si (10355) 


0.64 


0.44 


0.476 


0.263 


0.443 ± 0.008 


4'' Si (10579) 


0.51 


0.35 


0.248 


0.104 


0.272 ± 0.029 


5^ Si (10865) 


0.42 


0.29 


0.310 


0.04 


0.31 ±0.07 


6^Si(11019) 


0.37 


0.25 






0.130 ±0.030 


7''Si(11193) 


0.32 


0.22 









TABLE III: Two-photon decay widths (in units of eV) of the pseudoscalar (^So), scalar (^Po), and tensor (^P2) bottomonium 
states. Bottomonium masses are in units of MeV. 







Theory 


Experiment 


state 


mass 


Ref.124] 


Ref.f25] 


Rcf.f4] 


Ref.^ 


Ref.f27] 


Ref.[28l 


Ref.f29] 


Ours 


PDG [18] 




9389 


350 


220 


214 


266 


192 


460 


460 


527 






9987 


150 


110 


121 


95.0 


116 




200 


263 






10330 


100 


84 


90.6 


67.9 


93.5 






172 






10595 




71 


75.5 


56.3 


81.8 






105 






10817 


121 




Vb{6'So) 


11011 


50 






9859 


38 


24 


20.8 


27.3 


24.1 


80 


43 


37 




Xbo{2'Po) 


10233 


29 


26 


22.7 


26.9 


27.3 






37 




XboiS'Po) 


10502 


35 




Xb2{l'P2) 


9912 


8 


5.6 


5.14 


2.56 


6.45 


8 


7.4 


6.6 




Xb2{2'P2) 


10269 


6 


6.8 


6.21 


6.11 


8.1 






6.7 




Xb2{S'P2) 


10540 


6.4 
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TABLE IV: El transition rates of bottomonium states in our screened potential model (those calculated by the zeroth-order 
wave functions are marked by SNRq and those by the first-order relativistically corrected wave functions are marked by 
SNRi). We also list the results of one potential model, in which the potential is determined by the inverse-scattering method, 
for comparison [sol ]. 
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81 


92 


0.86 


1.26 


1.13 
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state 


Initial meson 


Final meson 




(MeV) 




Pthy (keV) 




Pexpt (keV) 








Ref [30] 


SNRoa) 


Ref [30] 


SNRo 


SNRi 


PDG [18] 


3P ^ ID 


X62(3^P2) 


T(1^D3)(10156) 




377 


^ 


fa 


0.05 








T(1^D2)(10151) 




381 





^ 


« 








T(l^Di)(10145) 




387 


^ 


« 


« 






Xh(3^Pi) 


T(1^D2) 




366 





r; 


0.09 












372 


^ 


fa 


0.004 






Xfco(3'Po) 


T(l^Di) 




351 


^ 


^ 


0.17 








/ii,2(l^D2) (10152) 




370 




^ 


0.24 




ID IP 


T(1^D3)(10156) 


X62(1'P2) 


245 


240 


24.3 


26.4 


24.5 






T(1^D2)(10151) 


Xf>2(l^P2) 


240 


236 


5.7 


6.29 


5.87 








Xm(I^Pi) 


261 


255 


22.0 


23.8 


19.8 






TDi)(10145) 


Xf>2(l^P2) 


233 


230 


0.58 


0.65 


0.61 








Xm(I'Pi) 


254 


249 


11.3 


12.3 


10.3 








X6o(l^Po) 


285 


282 


21.4 


23.6 


16.7 






/l62(l'D2)(10152) 


/l6(l'Pl) 




246 




42.3 


36.5 




2D 2P 


T(2^D3)(10442) 


X62(2^P2) 


174 


172 


16.3 


18.0 


15.9 






T(2^D2)(10438) 


Xf>2(2^P2) 


171 


168 


3.83 


4.17 


3.82 








Xm(2^Pi) 


183 


181 


14.2 


15.7 


12.1 






TDi)(10432) 


Xf>2(2^P2) 


165 


162 


0.38 


0.42 


0.39 








Xm(2^Pi) 


178 


175 


7.2 


7.87 


6.35 








X6o(2^Po) 


202 


198 


14.2 


15.1 


9.49 






/ii,2(2^D2)(10439) 


/i6(2iPi) 




181 




31.3 


25.4 




2D IP 


T(1^D3) 


X62(1'P2) 


518 


517 


3.94 


4.01 


3.73 








, , /I 3-n \ 
XMl^i r2) 


014 


olo 


n 07 
u.y / 


U.yo 


U.Do 








X6l(l^Pl) 


534 


531 


3.25 


3.26 


4.46 






TDi) 


Xi>2(l^P2) 


509 


507 


0.10 


0.11 


0.05 








X6l(l^Pl) 


529 


525 


1.75 


1.76 


1.87 








Xi>o(l^Po) 


559 


557 


2.76 


2.79 


6.20 






/l62(l'D2) 


/l6(l'Pl) 




522 




6.19 


7.30 





TABLE V: Hyperfine and fine splittings in units of MeV for charmonium and bottomonium in our model. Here a is 1.362 GeV for 
charmonium and 3.3 GeV for bottomonium. The experimental values are the mass differences of the corresponding charmonium 
and bottomonium states taken from PDG [l3]. Results of Ref. [30] and Ref. [4j are listed for comparison. 



State 




Charmonium 




Bottomonium 






Ours 


Exp 


Ours 


Ref. [36] 


Ref [4] 


Exp 


l^Si - I'So 


118 


116.6 ± 1.2 


71 


87 


60 


71.4+3 ? ±3.7 


2^Si - 2^5*0 


50 


52 ±4 


29 


44 


20 




3^Si -3'5'o 


31 


21 


41 


10 




l''P2 - l^Pi 


44 


45.54 ±0.11 


21 


22 


20 


19.43 ±0.57 


l^Pl - l^Po 


77 


95.91 ±0.32 


32 


30 


30 


33.34 ± 0.66 


2^P2 - 2^Pi 


36 


18 


18 


10 


13.19 ±0.77 


2-'Pi - 2^Po 


59 


25 


25 


20 


22.96 ±0.84 


3^P2 - 3^Pi 


30 


16 


3^ Pi - 3^Po 


47 


22 



